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Summary 
Major  Results 

(i ) A breakthrough  in  substorm  prediction 

It  is  of  vital  importance  to  the  USAF  to  predict  accurately  the 
occurrence  of  magnetospheric  substorms.  The  success  of  this  task  would 
depend  on  whether  or  not  one  could  find  an  interplanetary  quantity  which 
is  physically  meaningful  and  is  closely  related  to  the  substorm  index 
AE  It  has  already  been  known  that  the  north-south  component  (namely, 
the  component  perpendicular  to  the  ecliptic  plane)  of  the  interplanetary 
magnetic  field  plays  an  important  role  in  the  occurrence  of  substorms. 
However,  it  is  too  crude  to  be  a parameter  in  predicting  both  the  occurrence 
and  intensity  of  substorms.  For  this  reason,  we  have  made  an  intensive 
search  for  physically  meaningful  parameters. 

Knowing  that  the  daily  sum  of  the  Kp  index  (denoted  by  IKp)  is  well 
correlated  with  the  solar  wind  speed  (Snyder,  Neugebauer  and  Rao,  1963; 
Olbert,  1968)  and  that  the  Bz  component  has  an  important  role  in  triggering 
substorms  (Arnoldy,  1971)  and  that  the  variance  of  the  IMF  correlates 
with  the  Kp  index  (Ballif,  Jones  and  Col  man,  1969),  we  examined  also  the 
relationship  (i)  between  the  variance  of  Bz  and  the  AE  index  and  (ii) 
the  variance  of  Bz  and  solar  wind  speed.  However,  it  is  concluded  that 
one  can  hardly  say  that  a high  speed  variance  of  the  IMF  Bz  causes  an 
intense  substorm  activity.  This  statement  can  also  be  easily  demonstrated 
by  examining  the  relationship  between  the  variance  of  the  IMF  Bz  component 
and  the  IKp  index  as  a function  of  time. 


After  an  extensive  search,  however,  it  is  found  that  the  parameter 
e(t)  = VB2  (£q  sin2  f^)2 

is  most  closely  related  to  the  AE  index.  This  quantity  e(t)  is  the 

* 

amount  of  Poynting  flux  which  flows  into  the  magnetosphere,  where 

JdS  = S - {£0(1  - sin  0/2) } 2 
= {l0  sin2  |-}2 
iQ  = 7 earth  radi i , 

and  0 is  a measure  of  the  angle  between  the  geomagnetic  field  vector 
and  the  IMF  vector  at  the  front  of  the  magnetosphere  in  the  equatorial 
plane.  The  above  expression  for  _fdS  was  obtained  empirically  as  a 
result  of  an  effort  to  find  the  most  suitable  expression  which  relates 
the  Poynting  flux  to  the  rate  of  the  total  energy  dissipation  u(t). 

This  can  be  seen  in  the  relationship  between  e(t)  and  the  AE  index. 
Figure  la  shows  this  relationship  for  one  of  the  typical  recurrent 
storm:,  November  3_10,  1973.  This  figure  may  be  contrasted  with  Figure 
lb  which  shows  the  relationship  between  the  IMF  Bz  component  and  the  AE 
index  for  the  same  period;  there  is  hardly  a quant i tative  relationship 
between  them.  Figures  2a  and  2b,  and  3a  and  3b  show  similar  examples; 
for  more  examples,  see  the  attached  paper.  Recognizing  the  the  AE  index 
does  not  necessarily  represent  accurately  the  substorm  intensity,  the 
• obtained  correlation  is  very  assuring. 
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Fig  lb 
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It  is  important  to  note  in  the  expression  of  e that  the  size  of  the 
cross-sectional  area  is  controlled  by  the  relative  angle  9'  between  the 
interplanetary  and  magnetospheric  magnetic  field  vectors.  It  becomes 
null  only  when  9'  = 0.  In  the  past,  it  has  been  emphasized  that  a 
substorm  occurs  when  the  north-south  component  of  the  IMF  becomes  negative 
or  only  when  it  is  directed  southward.  It  is  our  finding  that  the  AE 
index  does  not  become  very  small  unless  9'  = 0.  Thus,  a relatively 
intense  substorm  can  occur  even  when  9'  = 45°  and  thus  the  north-south 
component  of  the  interplanetary  magnetic  field  is  even  positive. 

Furthermore,  in  Arnoldy's  empirical  relation, 

AE  = -0.26  (ZBst)0  - 0.91(IBsx)1  - 0.33  (ZBst)2  + 0.12  PQ, 

the  quantity  Bst  has  a value  only  when  Bs  = Bz  < 0 and  thus  the  last 
term  (practically  a constant)  is  introduced  to  make  a continuous  function. 
Further,  the  physical  meaning  of  the  first  three  terms  is  not  necessarily 
very  clear.  On  the  other  hand,  our  relation  represents  a definite 
physical  meaning,  namely  the  exact  amount  of  the  interplanetary  energy 
which  is  transported  into  the  magnetosphere  across  the  magnetopause 
which  is  subsequently  converted  into  the  substorm  energy.  Thus,  by 
monitoring  e(t)  in  interplanetary  space  at  the  lagrangian  point,  it  is 
possible  to  predict  both  the  occurrence  and  intensity  of  magnetospheric 
substorms  about  30-60  minutes  prior  to  their  occurrence. 
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We  feel  that  this  finding  is  a major  breakthrough  in  the  search  for 
physically  meaningful  parameters  which  can  quantitatively  predict  the 
occurrence  and  intensity  of  substorms. 

The  quantity  e(t)  is  also  tested  for  major  geomagnetic  storms  with 
a large  main  phase  decrease.  It  is  our  finding  that  the  quantity  e(t) 
agrees  over  three  orders  of  magnitude  with  the  total  amount  of  energies 
produced  within  the  magnetosphere,  confirming  its  accuracy  even  during 
violent  geomagnetic  activity.  Figures  4a-c  show  an  example  of  the 
comparison  between  t and  the  energy  dissipated  in  the  magnetosphere 
(Perreault  and  Akasofu,  1978). 

(ii)  Study  of  the  simultaneous  auroral  photographs  taken  from  above 
by  the  OMSP  satellite  and  from  below  by  the  all-sky  camera 
We  made  an  extensive  study  of  the  simultaneous  auroral  photographs 
taken  from  above  by  the  DMSP  satellite  and  from  below  by  the  all-sky 
camera  at  the  South  Pole  Station.  Figures  5a  through  5g  show  a complete 
set  of  data  for  June  12,  1975.  Each  figure  shows  a OMSP  photograph  and 
the  simultaneous  all-sky  photograph  taken  from  the  South  Pole.  The 
location  of  the  South  Pole  in  the  DMSP  photograph  is  indicated  by  a 
white  dot;  the  location  of  the  invariant  pole  is  indicated  by  a white 
cross;  the  noon  meridian  by  a white  line.  In  the  all-sky  photograph, 
the  direction  of  the  invariant  pole  is  indicated  by  a white  bar  near  the 
edge  of  the  circular  photograph.  Each  figure  shows  several  all-sky 
photographs  taken  shortly  before  and  after  the  satellite  passage.  Both 
OMSP  and  all -sky  photographs  complement  each  other  (the  former  provide  a 
large-scale  view,  and  the  latter  provide  fine-time-scale  data).  For 
more  details,  see  Akasofu  (1978). 
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Fig  5e 
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Other  Results 

(i ) Equatorward  boundary  of  the  electron  precipitation 

Electron  spectrograms  from  351  passes  of  the  Isis  1 and  2 satellites 
were  utilized  to  study  statistically  the  effects  of  the  interplanetary 
magnetic  field  (IMF),  substorm  activity,  and  the  earth's  dipole  tilt 
angle  on  the  latitude  of  the  equatorward  boundary  of  the  nightside 
(2000-0400  magnetic  local  time)  auroral  oval.  The  boundary  location  (in 
invariant  latitude)  at  hourly  local  time  intervals  was  identified  in 
terms  of  the  equatorward  boundary  of  the  diffuse,  > 100  eV  electron 
precipitation.  The  following  characteristics  were  noted:  (1)  The  north- 
south  component  (B^)  of  the  IMF  plays  the  dominant  role  in  controlling 
the  motion  of  this  boundary.  The  invariant  latitude  of  the  boundary  is 
shown  to  shift  by  approximately  + 4°  depending  on  the  direction  of  the 
IMF  (northward  and  southward,  respectively)  relative  to  its  position 
corresponding  to  Bz  = 0.  This  indicates  an  inward  motion  of  the  associated 
boundary  in  the  magnetotail  by  about  5 earth  radii  when  the  IMF  changes 
its  direction  from  northward  to  southward  with  large  magnitude.  There 
is  a significant  difference  in  the  amount  of  the  shift  between  the 
evening  and  morning  sectors,  i.e.,  for  the  same  decrease  in  82  value, 
the  boundary  moves  more  equatorward  in  the  morning  sector  than  in  the 
evening  sector.  When  the  obtained  oval  particle  boundary  is  projected 
onto  the  equatorial  plane  of  the  magnetotail  along  magnetic  field  lines, 
a good  agreement  was  found  between  the  projected  boundary  and  the  drift 
boundary  (the  Alfven  layer)  of  low  energy  electrons  in  the  presence  of 
the  dawn-dusk  electric  field.  Thus,  this  agreement  gives  new  evidence 
showing  that  the  diffuse  auroral  particles  originate  near  and  at  the 
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inner  boundary  of  the  plasma  sheet.  (2)  Substorm  activity  seems  to  have 
a separate  role  in  determining  the  latitude  of  the  equatorward  boundary 
of  the  nightside  auroral  precipitation  region.  The  boundary  during 
substorm  periods  is  statistically  found  to  be  2 - 3°  lower  in  invariant 
latitude  than  that  during  quiet  times.  Even  a simple  classification 
into  quiet  and  disturbed  conditions  improves  the  accuracy  with  which  the 
auroral  oval  location  can  be  inferred.  By  combining  the  IMF  effect  and 
the  substorm  effect,  it  is  indicated  that  the  boundary  is  located  at  the 
lowest  latitudes  when  a substorm  takes  place  during  a southward  IMF  with 
large  magnitude,  whereas  it  is  located  in  the  highest  latitudes  when  the 
IMF  has  a northward  component  during  quiet  times.  (3)  The  equatorward 
boundary  of  the  nightside  auroral  oval  is  located  in  higher  latitudes  in 
the  winter  hemisphere  than  in  the  summer  hemisphere,  although  this 
earth's  dipole  tilt  effect  is  usually  smaller  than  the  effects  of  the 
IMF  and  substorm  activity.  For  more  details,  see  Kamide  et  al.  (1978). 

( i i ) Correlation  study  between  the  power  of  the  meter-wave  radiation 
and  the  AE  index 

In  the  past  several  years  low  frequency  radio  measurements  have 

shown  that  the  earth's  magnetosphere  is  a very  intense  radio  emitter, 

with  characteristics  very  similar  to  other  astronomical  radio  sources 

such  as  Jupiter  and  Saturn.  Gurnett  showed  that  these  radio  emissions 

are  closely  associated  with  the  occurrence  of  auroral  arcs  on  the  night 

side  of  the  earth.  Both  the  angular  distribution  of  the  emitted  radiation 

and  direction  finding  measurements  indicate  that  such  an  intense  (up  to 

g 

10  watts,  total  radiated  power)  radiation  associated  with  auroral  arcs 
is  generated  at  relatively  low  altitudes  (^  3 radial  distance)  over 
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the  night  and  evening  auroral  regions,  presumably  by  the  intense  fluxes 
of  several  keV  electrons,  causing  the  aurora  and  the  currents  responsible 
for  the  magnetic  disturbances.  Because  of  the  close  association  with 
auroral  phenomena  and  the  wavelengths  in  the  kilometer  range  these  radio 
emissions  have  been  called  auroral  kilometric  radiation. 

Although  kilometric  radiation  is  known  to  be  associated  with  high 
latitude  magnetic  disturbances  little  has  been  done  to  study  this  relation- 
ship in  detail.  It  is  of  interest  to  determine  just  how  good  the  correlati 
is  and  whether  there  are  any  exceptions  to  the  observed  relationship. 

The  relationship  of  the  radio  emission  intensity  to  the  currents  flowing 
through  the  auroral  zone  may,  for  example,  be  helpful  in  developing  a 
better  understanding  of  how  the  auroral  kilometric  radiation  is  generated. 
Furthermore,  since  radio  emissions  from  the  entire  auroral  zone  can  be 
easily  monitored  by  a single  spacecraft  there  is  the  question  of  whether 
the  auroral  kilometric  radiation  intensity  could  provide,  on  a near  real 
time  basis,  a parameter  which  is  comparable  to  magnetic  disturbance 
indices,  such  as  the  auroral  electrojet  index,  AE,  which  require  a large 
array  of  ground  stations  with  the  attendant  problems  of  information 
retrieval  and  processing.  We  made  a detail  study  of  the  relationship 
between  the  auroral  kilometric  radiation  intensities  observed  by  the 
IMP-6  spacecraft  and  the  AE  index  computed  from  a series  of  high- 
latitude  magnetometer  station  and  demonstrated  that  the  auroral  kilometric 
radiation  is  a good  indicator  of  auroral  magnetic  disturbances  (Figures 
6a  and  6b);  for  more  details,  see  Voots  etal.  (1977). 
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(iii)  Auroral  circle 


The  accurate  determination  of  the  shape  of  the  auroral  oval  is  an 
important  problem  in  magnetospheric  physics,  as  well  as  in  predicting 
the  location  of  the  oval.  C.-I.  Meng,  R.H.  Holzworth  and  the  Principal 
Investigator  made  an  extensive  study  on  this  subject  by  analyzing  more 
than  fifty  DMSP  auroral  photographs.  It  is  found  that  the  poleward 
boundary  of  the  auroral  oval  can  be  fitted  by  an  off-centered  circle 
within  an  accuracy  of  less  than  0.5°  in  latitude;  for  more  details,  see 
Meng  et  al.  (1977). 

Civ)  Polar  cap  arcs 

: 

It  has  been  well  established  that  auroras  tend  to  appear  along  an 
oval  band  known  as  the  auroral  oval.  However,  there  is  another  type  of 
aurora  which  appears  well  inside  the  area  bounded  by  the  auroral  oval, 
namely  polar  cap  auroras.  C.-I.  Meng  and  the  Principal  Investigator 

i 

< 

have  found  several  DMSP  auroral  photographs  which  show  polar  cap  arcs 
which  lie  parallel  to  oval  arcs  by  changing  the  orientation,  from  the 
sun-earth  line  (or  the  noon-midnight  meridian)  direction  to  the  westward 
direction,  near  the  poleward  boundary  of  the  oval  in  the  night  sector; 
for  more  details,  see  Meng  and  Akasofu  (1976). 

(v)  Response  of  the  dayside  aurora  to  sharp  northward  and  southward 
transitions  of  the  interplanetary  magnetic  field  and  to 

- 

magnetospheric  substorms  ' 

Latitudinal  shifts  of  the  dayside  aurora  are  examined  with  respect 

1 

to  variations  in  the  interplanetary  magnetic  field  (IMF)  and  to  magneto- 
spheric substorms.  Within  10-15  minutes  after  steplike  southward  (northward) 
transitions  in  the  IMF,  the  dayside  auroral  oval  moves  equatorward 
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(poleward).  In  some  cases,  the  auroral  shift  is  similar  to  an  exponential 
relaxation  in  latitude  from  an  initial  to  a final  steady  state  value; 
for  these  cases,  the  average  exponential  time  constant  is  estimated  to 
be  17  minutes.  Substorm  features  in  the  dayside  aurora  include:  (a)  an 
equatorward  shift  by  1°  to  3°  at  the  equatorward  auroral  boundary, 

(b)  a brightening  of  the  aurora  near  substorm  onset,  (c)  the  formation 
of  multiple  auroral  bands  and  (d)  poleward  motion  of  individual  auroral 
arcs  occurring  nearly  coincident  in  time  with  the  equatorward  boundary 
shifts,  the  average  poleward  velocity  of  these  arcs  being  approximately 
800  m/sec. 

(vi)  HF  radio  wave  absorption 

It  has  been  well  known  that  auroral  activity  causes  a heavy  absorption 
of  HF  radio  waves.  However,  because  of  a lack  of  an  appropriate  array 
of  riometer  stations,  the  latitudinal  development  of  cosmic  noise  absorption 
during  auroral  substorms  has  not  been  accurately  known. 

The  operation  of  both  all-sky  camera  and  riometers  along  the  Alaskan 
meridian  in  1973  made  such  a study  possible.  Changes  of  the  latitudinal 
distribution  of  auroras  were  carefully  examined  over  a latitudinal  range 
of  more  than  10^,  and  the  simultaneous  changes  were  obtained  by  using  a 
computer  contouring  technique.  A paper  on  this  subject  is  now  near 
completion  (Berkey,  Anger  and  Akasofu). 

(vii)  Auroral  motions  and  the  (ExB)  drift  motion  of  ionospheric  plasma 

It  has  been  an  important  question  whether  or  not  auroral  motions 

represent  the  (ExB)  convective  motion.  Past  studies  suggested  that  the 
poleward  motion  of  auroras  during  the  expansive  phase  is  not  the  (ExB) 
motion,  but  the  eastward  motions  in  morning  hours  are.  However,  there 
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had  been  no  extensive  study  of  this  important  subject,  since  a systematic 
monitoring  of  the  ionospheric  electric  field  was  not  available  until 
1970,  when  the  Chatanika  incoherent  scatter  radar  came  into  existence. 

However,  during  the  last  several  years,  a large  amount  of  high 
quality  electric  field  measurements  and  the  simultaneous  all-sky  photographs 
became  available.  Therefore,  we  made  an  extensive  study  of  the  correlation 
between  the  two  data  sets.  The  results  are  shown  in  Figure  7.  Our 
main  result  is  that  auroral  motions  are  generally  far  from  the  (ExB) 
drift  motion  in  the  ionospheric  level.  Often,  they  are  opposite  to  the 
(ExB)  drift  motion. 
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Figure  la  Correlation  between  the  auroral  electrojet  index  AE  and 


the  interplanetary  quantity  e(t)  for  the  November  3-10, 
1973  storm. 

Figure  lb  Correlation  between  the  auroral  electrojet  index  AE 
and  the  Bz  component  of  the  interplanetary  magnetic 
field  (IMF)  for  the  November  3-10,  1973  storm.  The 
8Z  component  is  plotted  in  reverse  and  negative  periods 
are  hatched;  positive  periods  are  plotted  by  a dotted  line. 

Figure  2a  Correlation  between  the  AE  index  and  e(t)  for  the  April 
16-23  storm. 

Figure  2b  Correlation  between  the  AE  index  and  the  Bz  component 

of  the  IMF  for  the  April  16-23,  1974  storm;  for  details, 
see  the  caption  for  Figure  la. 

Figure  3a  Correlation  between  the  AE  index  and  e(t)  for  the 
January  22-29,  1974  storm. 

Figure  3b  Correlation  between  the  AE  index  and  the  Bz  component 
of  the  IMF  for  the  January  22-29  storm;  for  details, 
see  the  caption  for  Figure  la. 
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Figure  4a  The  magnetic  field  !BI  and  Bz  and  the  combined  auroral 

zone  magnetic  record  (AE*)  and  the  combined  low  latitude 
magnetic  record  (ASY*)  for  the  storm  of  February  7-9,  1967. 

Figure  4b  Comparison  of  the  interplanetary  quantities  |B|,  Bz, 

ZBjt,  V x Bz  (westward  component  only)  and  the  magnetospheric 
quantities  AE,  Dst,  ASY,  AL  for  the  storm  of  February 
7-9,  1967. 

Figure  4c  Comparison  of  the  interplanetary  quantities  e and  the 
energy  dissipated  in  the  magnetosphere  which  includes 
both  the  joule  heat  dissipation  and  the  ring  current 
injection. 

I 1 

Figures  5a-g  DMSP  photograph  of  the  antarctic  oval  on  June  12,  1975 
and  the  simultaneous  all -sky  photograph  taken  from  the 
South  Pole:  A dot  with  an  arrow  in  the  DMSP  photograph 
indicates  the  location  of  South  Pole.  In  the  all-sky 
photograph,  a white  bar  indicates  an  approximate  location 
of  the  invariant  pole.  In  the  lower  half  of  the  figure, 
selected  all-sky  photographs  from  South  Pole  are  shown, 
indicating  changes  of  the  auroras  at  about  the  time  of 
the  satellite  passage. 


Figure  6a  A scatter  plot  of  simultaneous  measurements  of  the 

auroral  kilometric  radiation  intensity  (normalized  to 
30  R^)  and  the  auroral  electrojet  index  AE  using  a 
random  selection  of  50  days  of  data  from  the  IMP-6 
spacecraft. 

Figure  6b  Simultaneous  measurements  of  the  auroral  kilometric 
radiation  intensity  at  178  kHz  by  the  IMP-6  satellite 
and  the  auroral  electrojet  index,  AE,  for  a selection  of 
days  in  which  a close  correlation  is  observed.  The 
coordinates  of  IMP-6  are  shown  at  the  bottom  of  each 
panel. 

Figure  7 Comparison  of  the  observed  north-south  component 

electric  field  by  the  Chatanika  radar  and  the  inferred 
north-south  electric  field  from  the  east-west  auroral 
motions.  If  auroral  motions  would  arise  from  the 
ionospheric  electric  field,  all  points  should  align  along 
the  line  bisecting  the  first  and  fourth  quadrants. 
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Abstract 

In  Part  II,  our  new  finding  of  the  relationship  between  e(t)  and 
the  AE  index  will  be  described  in  more  detail.  In  particular,  it  is 
shown  that  the  interplanetary  quantity  e(t),  obtained  by  Perreault  and 
Akasofu  (1978),  for  intense  geomagnetic  storms,  also  correlates  well 
with  individual  magnetospheric  substorms.  This  quantity  is  given  by 
e ( t ) * V | B | ^ sln^  (e/2)  tQ2,  where  V and  |B[  denote  the  solar  wind  speed 
and  the  magnitude  of  the  interplanetary  magnetic  field  (IMF),  respectively, 
and  9 denotes  the  polar  angle  of  the  IMF;  is  a constant  * 7 earth 
radii.  The  AE  index  is  used  in  this  correlation  study.  The  correlation 
is  good  enough  to  predict  both  the  occurrence  and  intensity  of  magnetospheric 
substorms  observed  in  the  auroral  zone,  by  monitoring  the  quantity  e(t) 
upstream  of  the  solar  wind. 
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1.  Introduction 


A number  of  studies  have  been  made  in  the  past  in  an  attempt  to 
find  interplanetary  quantities  which  are  related  to  magnetospheric 
substorm  activity.  Snyder,  Neugebauer  and  Rao  (1963)  and  Olbert  (1968) 
found  a good  correlation  between  the  daily  sum  of  the  Kp  index  (denoted 
by  ZKp)  and  the  solar  wind  speed.  It  is  given  by 

V( km/sec)  =6.3  ZKp  + 262. 

Ballif,  Jones  and  Coleman  (1969)  showed  that  the  variance  of  the  inter- 
planetary magnetic  field  is  correlated  with  the  Kp  index.  Fairfield  and 
Cahill  (1966)  and  Rostoker  and  F^lthammar  (1967)  were  the  first  to 
recognize  the  importance  of  the  Bz  component  of  the  interplanetary 
magnetic  field  (IMF)  on  the  accurrence  of  magnetospheric  substorms. 

Since  then  a large  number  of  papers  have  been  published,  reporting  a 
good  correlation  between  the  southward  component  of  Bz  (or  Bz  < 0)  and 
the  magnetospheric  substorm  index  AE  (cf.  Arnoldy,  1971;  Burton  et  al , 

1975:  for  reviews,  see  Burch,  J.L. , 1974;  Akasofu,  1977). 

Most  recently,  however,  Perreault  (1974),  Burlaga  and  lepping 
(1977)  and  Perreault  and  Akasofu  (1978)  have  examined  closely  this 
correlation  and  have  found  that  the  dependence  of  substorm  activity  on 
the  Bz  component  is  far  from  satisfactory.  Burlaga  and  Lepping  (1977) 
pointed  out  the  importance  of  the  solar  wind  speed  in  the  correlation, 
in  particular  in  a later  period  of  recurrent  geomagnetic  storms.  Perreault 
(1974)  and  Perreault  and  Akasofu  (1978)  noted  that  a significant  substorm 
activity  can  occur  even  when  the  Bz  component  is  positive  and  thus  that 
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the  magnetosphere  cannot  be  a rectifier  which  passes  only  negative 


values  of  the  Bz  component.  They  found  also  that  the  magnitude  of  the 
Bz  component  alone  cannot  predict  quantitatively  the  intensity  of  magneto- 
spheric  substorms  and  storms. 

In  their  search  for  interplanetary  quantities  which  can  provide  a 
better,  quantitative  correlation  with  substorm  activity,  they  found  the 
quantity  e(t)  given  by 

e(t)  = V|B|  2 sin4  f lZ 

where 

V = the  solar  wind  speed 

|Bt  = the  magnitude  of  the  IMF 

0 = the  polar  angle  of  the  IMF  vector  in  the  y-z  plane 

in  solar-magnetospheric  coordinates 

= tan'1  (|  B | /Bz)  for  Bz  > 0 

= 7t  - tan  1 (|By|/Bz)  for  Bz  < 0 

2,  = 7 earth  radii, 
o 

The  correlation  between  e(t)  and  the  total  energy  U(t)  generated  in  the 
magnetosphere  during  twelve  major  geomagnetic  storms  is  found  to  be 
quite  good.  The  total  dissipated  energy  U includes  the  ring  current 
energy,  joule  heat  energy  associated  with  the  auroral  electro jet  and  the 
kinetic  energy  deposited  in  the  polar  ionosphere  by  auroral  particles. 
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However,  since  all  the  major  geomagnetic  storms  represent  periods 
during  which  intense  substorms  occur  very  frequently,  it  was  not  possible 
for  Perreault  and  Akasofu  (1978)  to  examine  in  detail  the  relationship 
• between  s(t)  and  individual  substorms. 

The  purpose  of  this  paper  is  to  examine  the  relationship  between 
e(t)  and  the  AE  index  for  weak  storms  during  which  individual  substorms 
can  reasonably  be  resolved.  Fortunately,  the  necessary  interplanetary 
quantities  in  estimating  e(t)  are  available  for  some  of  the  recurrent 
geomagnetic  storms  in  the  interval  1973-1975.  All  quantities  dealt  with 
in  this  paper  are  hourly  average  values. 

2.  Results 

(a)  The  November  3-10.  1973  Storm 

This  particular  storm  was  examined  by  Burlaga  and  Lepping  who 
noted  that  a good  quantitative  correlation  between  the  Bz  component  and 
the  AE  index  holds  only  during  an  initial  period  of  recurrent  storms  and 
that  the  solar  wind  speed  plays  an  important  role  in  a later  period. 

In  Figure  la,  the  relationship  between  the  AE  index  and  the  8Z 
component  is  shown  for  this  particular  storm.  The  Bz  component  was 
plotted  in  reverse;  a negative  period  is  hatched  and  a positive  period 
is  plotted  by  a dotted  line.  Further,  the  two  quantities  are  normalized 
in  such  a way  that  peaks  of  both  quantities  have  similar  'heights' 
during  an  initial  period  of  the  storm.  As  recognized  correctly  by 
Burlaga  and  Lepping  (1977),  a good  correlation  breaks  down  after  the 
second  day  of  the  recurrent  storm. 
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In  Figure  lb,  e(t)  and  the  AE  index  are  plotted  by  a solid  line  and 

a dotted  line,  respectively.  One  can  see  immediately  considerably 

higher  correlation  between  e(t)  and  the  AE  index  and  between  B,  and  the 

z • 

AE  index.  Note  in  particular  that  e(t)  traces  out  reasonably  well  AE 
index  variations.  In  fact,  the  correlation  coefficient  for  the  e(t)-AE 
relation  is  0.665,  while  it  is  -0.477  for  the  BZ~AE  relation  for  this 
particular  storm.  For  reference,  the  correlation  coefficient  for  the 
VBZ“AE  correlation  is  0.459.  Since  the  AE  index  is  only  a rough 
index  for  substorm  activity  (Akasofu,  1977),  the  coefficient  value  of 
0.665  should  be  considered  to  be  quite  significant. 

(b)  The  June  7-14,  1974  Storm 

In  Figure  2a,  one  can  see  an  example  similar  to  (a). 

Figure  2a  also  illustrates  that  the  reasonably  good  correlation  between  Bz 
and  the  AE  index  breaks  down  after  the  middle  of  June  12  and  there  is 
little  quantitative  relationship  during  a later  period  of  the  storm.  On 
the  other  hand,  the  correlation  between  e(t)  and  the  AE  index  remains 
remarkably  high  throughout  the  period  (Figure  2b).  The  only  serious 
discrepancy  between  e(t)  and  the  AE  index  occurred  near  the  end  of  June 
10  and  the  beginning  of  June  11.  The  correlation  coefficients  between  e(t)  and 
AE,  Bz  and  AE,  and  VBZ  and  AE  are  0.639,  -0.580  and  -0.609,  respectively. 

I 

(c)  The  April  16-23.  1974  Storm 

This  is  another  good  example  to  show  the  breakdown  of  a good 
correlation  between  Bz  and  AE  during  a later  period  of  a storm  (Figure 
3a).  The  correlation  coefficient  between  Bz  and  AE  is  only  -0.352.  On 


the  other  hand,  the  correlation  between  e(t)  and  AE  is  0.565,  much 
higher  than  the  former  value.  For  reference,  the  correlation  coefficient 
between  VB2  and  AE  is  -0.447  for  this  particular  storm.  One  should  note 
also  that  in  this  particular  storm  the  correlation  between  e(t)  and  AE 
become  less  than  satisfactory  during  a later  period  of  the  storm. 

(d)  The  January  22-29,  1974  Storm 

This  particular  example  provides  an  additional  feature  which 
is  not  apparent  in  the  earlier  ones.  Again,  one  can  easily  recognize 
that  the  good  correlation  between  B2  and  AE  lasted  for  less  than  one  day 
and  there  is  little  quantitative  relation  during  the  rest  of  the  storm 
(Figure  4a).  The  correlation  coefficient  between  B2  and  AE  is  only 
-0.367.  The  correlation  between  e(t)  and  AE  is  significantly  higher, 
0.562,  than  the  former  value. 

However,  in  spite  of  a generally  good  correlation,  one  serious 
breakdown  occurred  on  January  25  during  which  e(t)  was  considerably 
greater  than  what  is  predicted  by  the  correlation  with  AE.  This  storm  is 
quite  unique,  compared  with  the  earlier  examples,  in  that  a significant 
main  phase  (-65  y)  developed  on  this  day.  It  appears  that  the  excess 
energy  flux  was  deposited  in  the  ring  current  belt.  We  shall  discuss 
this  feature  later. 

(e)  The  July  3-10.  1974  Storm 

This  storm  is  similar  to  the  previous  example  in  that  a large 
excess  of  e(t)  over  AE  occurred  during  the  storm.  Again,  the  correlation 
between  B2  and  AE  is  poor  except  for  the  first  three  days  (Figure  5a), 
while  the  correlation  between  c(t)  and  AE  is  reasonably  good  for  the 


entire  period.  The  only  exception  is  a large  excess  of  e(t)  during  an 
early  part  of  June  6.  In  fact,  an  intense  main  phase  of  -198  y developed 
during  the  same  period. 

3.  Discussion 

Perreault  (1974)  and  Perreault  and  Akasofu  (1978)  have  already 
shown  that  e(t)  correlates  well  with  the  total  energy  dissipation  U 
which  includes  the  ring  current  injection  and  the  kinetic  energy  flux  of 
auroral  particles,  as  well  as  the  joule  heat  production  by  the  auroral 
electrojet.  In  this  paper,  it  is  shown  that  e(t)  correlates  well  also 
with  the  AE  index  for  recurrent  storms  which  are  much  weaker  than  those 
examined  by  Perreault  and  Akasofu  (1978).  In  fact,  individual  substorms 
observed  at  auroral  zone  stations  are  well  correlated  with  peaks  of 
e(t). 

As  9.  first  approximation,  the  AE  index  is  proportional  to  the  rate 
of  jv  * ‘ . heat  energy  production  by  the  auroral  electrojet,  since  the 
variability  of  the  conductivity  of  the  ionosphere  is  much  greater  than 
that  of  the  electric  field.  The  'normalization'  employed  in  constructing 
the  e(t)-AE  correlation  figures  suggests  that  an  AE  value  of  1000  y corresponds 
to  an  input  value  e(t)  of  order  6.5  x 10*7  ergs/sec.  In  their  paper, 

Perreault  and  Akasofu  (1978)  obtained  the  joule  heat  production  of  5.5  x 
1017  ergs/sec  for  a 1000  y value  of  the  AE  index. 

The  quantity  V|B|  is  approximately  the  interplanetary  energy  flux 

2 2 • 

per  unit  area  and  ( lQ  sin  (8/2))  may  be  considered  to  be  the  cross-section 

through  which  the  energy  flux  flows  into  the  magnetosphere.  It  is  interesting  to 
note  that  V and  IB|  vary  much  more  smoothly  than  sin4  (0/2).  Thus,  the  term 
sin4  (6/2)  tends  to  construct  individual  substorms.  However,  all  the 
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three  parameters,  V,  |B|  and  sin4  (9/2),  contribute  in  determining  the 
intensity  of  substorms  observed  at  the  standard  auroral  zone  stations. 

It  is  also  important  to  recognize  that  the  correlation  with  the  AE  index 
. is  higher  with  e(t)  = V | B|  2 sin4  (0/2)  2Q2  than  with  B2  or  VBZ. 

Recognizing  that  the  AE  index  is  only  a rough  index  for  magnetospheric 
substorm  activity,  the  correlation  coefficient  between  e(t)  and  AE  of 
order  0.6  can  be  considered  to  be  satisfactory.  Indeed,  it  may  be 
worthwhile  to  re-examine  the  AE  index  in  some  instances.  In  particular, 
since  the  AE  index  cannot  properly  reflect  substorm  activity  along  a 
contracted  oval,  the  correlation  refers  only  to  substorms  observed  at 
the  standard  auroral  zone  stations. 

It  is  also  shown  that  when  e(t)  'exceeds'  greatly  the  energy  dissipated 
in  terms  of  the  joule  heat  production  in  the  ionosphere,  a significant 
ring  current  develops  in  the  magnetosphere.  Indeed,  it  has  been  shown 
by  a number  of  workers  that  the  rate  of  energy  injection  in  building  up 
an  appreciable  ring  current  exceeds  the  rate  of  joule  heat  energy  production 
(cf.  Akasofu,  1977,  p.  274). 

4.  Conclusion 

With  this  success  of  obtaining  a reasonable  correlation  between  e(t) 
and  AE,  we  should  be  able  to  predict  both  the  occurrence  and  intensity 
of  magnetospheric  substorms,  as  well  as  the  growth  of  the  ring  current, 

, if  the  interplanetary  quantities  V,  |B|,  and  6 can  be  monitored  on  a 

real  time  basis. 
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Figure  Captions 

Figure  la  Correlation  between  the  auroral  electrojet  index  AE 
and  the  Bz  component  of  the  interplanetary  magnetic 
field  (IMF)  for  the  November  3-10,  1973  storm.  The 
8Z  component  is  plotted  in  reverse  and  negative  periods 
are  hatched;  positive  periods  are  plotted  by  a dotted 
line. 

Figure  lb  Correlation  between  the  auroral  electro jet  index  AE  and 
the  interplanetary  quantity  e(t)  for  the  November  3-10, 
1973  storm. 


Figure  2a  Correlation  between  the  AE  index  and  the  Bz  component  of 
the  IMF  for  the  June  7-14,  1974  storm;  for  details  see 
the  caption  for  Figure  la. 


Figure  2b  Correlation  between  the  AE  index  and  e(t)  for  the 
June  7-14,  1974  storm. 

Figure  3a  Correlation  between  the  AE  index  and  the  Bz  component  of  the 
IMF  for  the  April  16-23,  1974  storm;  for  details,  see 
the  caption  for  Figure  la. 
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Figure  3b 


Figure  4a 


Figure  4b 


Figure  5a 


Figure  5b 


Correlation  between  the  AE  index  and  e(t)  for  the  April 
16-23  storm. 

Correlation  between  the  AE  index  and  the  Bz  component 
of  the  IMF  for  the  January  22-29  storm;  for  details 
see  the  caption  for  Figure  la. 

Correlation  between  the  AE  index  and  e(t)  for  the 
January  22-29,  1974  storm. 

Correlation  between  the  AE  index  and  the  Bz  component 
of  the  IMF  for  the  July  3-10,  1974  storm;  for  details, 
see  the  caption  for  Figure  la. 

Correlation  between  the  AE  index  and  e(t)  for  the 
July  3-10,  1974  storm. 
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